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We investigate the chiral phase transition at finite temperature (T) in colour SU(A^( = 3) Quan- 
tum Chromodynamics (QCD) with six species of fermions {Nj = 6) in the fundamental repre- 
sentation [|l|]. The simulations have been performed by using lattice QCD with improved stag- 
gered fermions. The critical couplings jSj^ for the chiral phase transition are observed for several 
temporal extensions A^,, and the two-loop asymptotic scaling of the dimensionless ratio Tc /Al 
(Al = Lattice Lambda-parameter) is found to be achieved for N > 6. Further, we collect jSjJ^ at 
Nf — (quenched), and Nf =4 at a fixed N =6 as well asNf = 8 at A'r =6,12, the latter relying 
on our earlier study. The results are consistent with enhanced fermionic screening at larger Nf. 
The ratio /Al depends very mildly on Nf in the Nf — — 4 region, begins increasing at Nf = 6, 
and significantly grows up at Nf — 8, as Nf reaches to the edge of the conformal window. We 
discuss the interrelation of the results with preconformal dynamics in the light of a functional 
renormalization group analysis. 
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1. Introduction 

Emergence of a conformal symmetry and a preconformal (walking) behavior in strongly fla- 
vored non-Abelian gauge theories has received much attention. Walking dynamics near the infra- 
red fixed point has been advocated as a basis for strongly interacting mechanisms of electroweak 
symmetry breaking. Lattice Monte-Carlo simulations are expected to provide a solid theoretical 
base to understand the (pre-)conformal nature in the gauge theory. 

A second zero of the two-loop beta-function of massless QCD with Nf flavours implies, at least 
perturbatively, the appearance of an infrared fixed point (IRFP) at A'y > 8.05 [Q] with the restoration 
of conformal symmetry before the loss of asymptotic freedom (LAP) at A^^^^ = 16.5. Conformality 
should emerge when the renormalized coupling at the would be IRFP is not strong enough to break 
chiral symmetry. This condition provides the lower bound Nj of a so called conformal window 
in the flavor space, and we find elaborated analytic predictions [^, ^ : for instance, the functional 
renormahzation group method ^ suggests N'^ ~ 12. Before the emergence of conformal symmetry, 
a qualitative change of dynamics is claimed ntNf = 6 based on instanton study [^. 

Recent lattice studies focused on the computation of the edge of the conformal window Nf 
and the analysis of the conformal window itself, either with fundamental fermions [^-|l^]. or other 



representations []16|]. Among the many interesting results with fundamental fermions, we single out 
the observation that QCD with three colours and eight flavours is still in the hadronic phase [Q, p^ , 
while Nf = 12 seems to be close to Nf'^, with some groups favouring conformality [||, ^ [iT], p^ ], 
and others chiral symmetry breaking [p^]. The onset of new strong dynamics at Nf = 6 has been 
implied via an enhancement of the ratio of chiral condensate to cubed pseudoscalar decay constant 
0- 

Using the thermal transition as a tool for investigating preconformal dynamics has been largely 
inspired by a renormahzation group analysis The critical temperature for the chiral phase 
transition has been obtained as a function of Nf. Then the onset of the conformal window has been 
estimated by locating the vanishing critical temperature. The phase transition line is almost linear 
with Nf for small Nf, and clearly elucidates the universal critical behaviour at zero and non-zero 
temperature in the vicinity of Nj-. Thus, it would be a promising direction to extend the knowledge 
of finite T lattice QCD to the larger Nf region, by using the FRG results as analytic guidance. 

In this proceedings, we investigate the thermal chiral phase transition for Nf = 6 colour 
SU(A'^c = 3) QCD by using lattice QCD Monte Carlo simulations with improved staggered fermions 
based on our recent study [[T]]. A^y = 6 is expected to be in the important regime as suggested by the 



results in Refs. [g, |8|]. We also compute the critical couplings for Nf = (quenched) and Nf = 4 at 



Nt = 6, and use the results from Ref. [10] for Nf = 8. Then we investigate Nf dependences of the 
chiral phase transition. 



2. Simulation setups 



Simulations have been performed in the same as in the study used for Nf = 8 in Ref. |10]: 



We have utilized the publicly available MILC code [18] with the use of an improved version of the 
staggered action, the Asqtad action, with a one-loop Symanzik [|19, 20] and tadpole [ p2[ ] improved 
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gauge action. The tadpole factor mq is determined by performing zero temperature simulations on 
the 12^ lattice, and used as an input for finite temperature simulations. 

To generate configurations with mass degenerate dynamical flavours, we have used the rational 



hybrid Monte Carlo algorithm (RHMC) [gJJ. Simulations for Nf = 6 have been performed by using 
two pseudo-fermions, and subsets of trajectories for the chiral condensates and Polyakov loop have 
been compared with those obtained by using three pseudo-fermions with the same Monte Carlo 
time step dT and total time length T of a single trajectory. We have observed very good agreement 
between the two cases for both evolution and thermalization. We have monitored the Metropolis 
acceptance and reject ratio, and adjusted T = 0.2 — 0.24 and dx = 0.008 — 0.016 to realize the best 
performance. 

Measured observables are the expectation values of the chiral condensate and Polyakov loop, 

A^f / r n\ 1 ' ^' 

J / T'..rA/r-il ^ 



where Ns (Nf) represents the number of lattice sites in the spatial (temporal) direction, 6^4 ,x is the 
temporal link variable, and tr^ denotes the trace in colour space. The output of this measurement is 
the critical coupling for the chiral phase transition. 

3. Results 

All results have been obtained for a fermion bare lattice mass am = 0.02. In the left panel of 
Figs. |l|, the expectation values of the chiral condensate a^{\lf\if) are displayed as a function of jSl 
for several Nf. It is found that different Nt give a different behaviour of a^{\ff\if)- The asymptotic 
scaling analysis below will confirm that it corresponds to a thermal chiral phase transition (or 
crossover) in the continuum limit. 

All values of the critical lattice coupling are summarized in Table |l|. For larger Nt, the 
signal for the chiral phase transition becomes less clear, hence we investigate the histogram of the 
chiral condensate: The histogram for A^f = 8 exhibits the double-peak structure at /3l = 5.2, i.e., the 
competition between chirally symmetric and broken vacua. The critical coupling can be estimated 
as ^ = 5.225(25) for N = 8. For N = 12, we also observe the double-peak structure in the 
histogram of the chiral condensate around jSl = 5.45. 

These results can be analyzed and interpreted in terms of the two-loop asymptotic scaling. Let 
us consider the two-loop lattice beta function, 

li{g) = -{bog'+hg') , (3.1) 
(^0, bi) = {{n-2Nf/3)/{4n)\ {m-3SNf/3)/{4nf) , (3.2) 

for fundamental fermions in colour SU(3). From Eq. ( |3.lD , we obtain the well known two-loop 
asymptotic scaling, 

Al a(/3L) = {2N,bo/PL) "''^^'""^ exp [-^l/(4A^c^o)] • (3.3) 

Here, Al is the so-called lattice Lambda-parameter, and pi = 2Nc/g^, with g = ^y2N^JlO ■ gi^. This 
definition effectively takes account of the improvement of the staggered lattice action when com- 



paring to the asymptotic scaling law, see Ref. |1C]. We insert Al to the definition of temperature 
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T=[ail5i^)Nt]-\ 

Nr' = (7;/Al) X (Al a(/3L^)) , (3.4) 
and extract the physical quantity Tc/Al by substituting the simulation outputs for Eq. (^). 



This ratio must be unique as long as the asymptotic scaling Eq. (3.3) is verified for a given 

In the right panel of Fig. |l|, the slope of the line connecting the origin and the data points corre- 
sponds to Tc/Al. The Nt = 6, 8, and 12 points have a common slope to a very good approximation, 
while the Nt = 4 result falls on a smaller slope. The latter is interpreted as a scaling violation effect 
due to the use of a too small A^,. The existence of a common T^/Al for Nt >6 indicates that the 
data are consistent with the two-loop asymptotic scaling Eq. ( |3.3| ), confirms the thermal nature of 
the transition and that Nf = 6is outside the conformal window, as expected from a previous Nf = S 
study A linear fit provides TjAi^ = 1.02(12) x 10^ which can be interpreted as the value in 
the continuum limit for Nf = 6 QCD. 

In order to have a more complete overview, we have performed simulations for the theory with 
Nf = (quenched) and Nf = 4, only at N = 6. These theories are of course very well investigated, 
however we have not found in the literature results for the same action as ours. We note that in a 



previous lattice study with improved staggered fermions []23|], asymptotic scaling was observed for 
A^, > 6 for < Nf < 4. Table | shows a summary of our results for the critical coupling fi^^ of the 
chiral phase transition at finite temperature for Nf = 0, 4, 6, and 8 - the latter from Ref. [|Io|]. 
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Figure 1: Left: The chiral condensate a^{\jf\i/) for Nf — 6 and am = 0.02 in lattice units, as a function of 
j3l, for A^, = 4, 6, 8, and 12. Error-bars are smaller than symbols. Right: The thermal scaling behaviour of 
the critical lattice coupling j3j^. Data points for Al a{Pi^) at a given 1/Nt are obtained by using (5^ from 
Table |l]as input for extracting Al a{l^) in the two-loop expression Eq. ( |3.3| ). The dashed line is a linear fit 
with zero intercept to the data with A', > 4. 

In the left panel of Fig. we display the critical values of the lattice coupling gc = ^JINcj^ 
from Table H in the Miransky-Yamawaki phase diagram. Consider the Nt = () results: it is expected 
that an increasing number of flavours favors chiral symmetry restoration. Indeed, we find that, on 
a fixed lattice, the critical coupling increases with Nf in agreement with early studies and naive 
reasoning. The precise dependence of the critical coupling on Nf at fixed Nt is not known. It is, 
however, amusing to note that the results seem to be smoothly connected by an almost straight line: 
the brown line in the plot is a linear fit to the data. Comparing the trend for Nf = 6 to the one for 



4 



Thermodynamic Study for Conformal Phase in Large Nf Gauge Theory 



Kohtaroh Miura 



Nf = ?i for varying Nt, one can infer a decreasing in magnitude (and small) step scaling function, 
hence a walking behaviour. Further study is needed at larger Nf, and by using the same action used 
for Nf = — S, to confirm or disprove it. 

Table 1: Summary of the critical lattice couplings ji^ for the theories with A^^ ~ 0, 4, 6, 8, am — 0.02 and 
varying A'r = 4, 6, 8, 12. All results are obtained using the same lattice action. 



Nf\N 
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7.88 ±0.05 
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4.675 ±0.025 
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FitN,=6 

N,=4 1 1 






















n1=8 I ■ 1 










N 


,= 12 
































































r4 



















0.85 0.9 0.95 



1.15 1.2 1.25 



3.5 



2.5 



1.5 



0.5 



B„^1.0 1-^— , 














( 


)- 












































J ^ 4 


t 



g,=(2VPL')"^ 



4 



Figure 2: Left: Critical values of the lattice coupling gc = y/2Nc/p^ for theories with Nf = 0, 4, 6, 8 
and for several values of A^, in the Miransky-Yamawaki phase diagram. The dashed (brown) line is a linear 
fit to the Nt = 6 results. Right; The Nf dependence of R{Nf)/R{0) for several finite fixed jSJ*^*. Here, 
R{Nf) = {Tc/Ar^f){Nf). 



Next, we study the Nf dependence of the ratio Tc/Al and related quantities. In addition to the 
scale Al, we introduce more UV reference energy scale A^ef, which is associated with a reference 
coupling jSl"^*. Then Eq. (33) is generalized as 



hi/{2bl) 



bl ^^-' + 2NM/bo 



exp 



^Ncbo 



ref 



(3.5) 



At leading order of perturbation theory b\ 0, we find A^f/AL = exp [jSj^'''^^/ (4A'^c^o)]- This equa- 
tion would be analogous of the ratio Al/Ams derived in [ p4| ] for Wilson fermions up to a further 
linear dependence on Nf in the numerator of the exponent. In a nutshell, the difference orig- 
inates from the fact that we are fixing a bare reference coupling jSJ*^^, which will be specified 
later. Notice that by construction Aref reproduces the lattice Lambda-parameter Al in the limit 

A,ef(i(3L'-'=*'^0)=AL(l±^(l/i8L=)). 
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Let us consider first R{Nf)\p ref^QQ = Tc/Al- The values of Tc/Al are found to be 600 ± 34, 
620 ±28, 1023 ±117, and 2098 ± 191 for Nf = 0, 4, 6, and 8, respectively, and represented as 
circles in the right panel of Fig. ^ (the vertical axis is normalized by R{0) = {Tc/Ai^){Nf = 0) for 
each Pi^^^)- The ratio does not show a significant Nf dependence in the region <Nf < 4, it starts 
increasing at Nf = 6, and undergoes a rapid rise around Nf = 8. The nearly constant nature of 
Tc/Al in the region Nf <4 indicates that the role of such energy scale is not significantly changed 



by the variation of Nf (see [ )25| ] for a detailed discussion of this point.) In turn, the increase of 
Tc/Al in the region Nf > 6 might well imply that the chiral dynamics becomes different from the 



one for Nf < 4. Indeed, a recent lattice study [ 17] indicates that Nf = 6 is close to the threshold for 
preconformal dynamics. 

We now consider Tc/Aref with finite The Nf dependence of the ratio R{Nf) = {Tc/ A^^f) (Nf) 
is shown for several jSj*^^ in the right panel of Fig. ^ (with normalization by R{0) = {Tc/ Aj:f.f){Nf = 
0) for each jSj*^^). Tc/Aj-^i is now a decreasing function of Nf for a larger fi^^^^. The Aref associ- 
ated with a pj^^^ ^ fi^ = 2Nc/g}^^ would be less sensitive to the IR or chiral dynamics. Assuming 
N'f ~ 12, the two-loop beta-function leads to jS* = —lNcb\/bQ ~ 0.63. The decreasing nature of 
(rc/Aref)(A^/) is found to start around ^'^'^ = 1.0 > j3*. Thus, the use of a UV reference scale 
leads to the decreasing {Tc/Arsf){Nf). This trend is consistent with the FRG study [^], where 
the decreasing Tc{Nf) has been obtained by using the T-lepton mass as a common UV ref- 
erence scale with a common coupling as{mj). We note that we have constrained our analyses 
jSj"^ < i3uv = Pi!{Nf) < 4.1 125 ± 0.0125. 

With the use of a UV reference scale, we should observe the predicted critical behavior [Q], 

Tc{Nf)=K\Nf-N}\-^l^ . (3.6) 

By choosing the critical exponent d in the range predicted by FRG: 1.1 < 1/|0| < 2.5, our data are 
consistent with the values A^j: = 9 ( 1 ) for ^^^^ = 4.0 and N'f =ll{2) for ^^"^ = 2. We plan to extend 
and refine this analysis in the future, and here we only notice a reasonable qualitative behaviour. 



4. Summary 

We have studied the chiral phase transition at finite T for colour SU (3) QCD with Nf = 6 
by using lattice QCD Monte-Carlo simulations with improved staggered fermions [|l]]. We have 
determined the critical lattice coupling ^ for several lattice temporal extensions Nt, and extracted 
the dimensionless ratio Tc/Al (Al =Lattice Lambda-parameter) by using two-loop asymptotic 



scaling. The analogous result for Nf = 8 has been extracted from Ref. [|10||. T^/Al for Nf = and 
Nf = 4 has been measured at fixed N = 6, barring asymptotic scaling violations. Then we have 
discussed the Nf dependence of the ratios Tc/Al and Tc/Aref, where Aref is a UV reference energy 
scale, related to Al via Aief/AL ~ exp[j3L'^'^^/(4A'^c^o)]- We have observed that Tc/Al shows an 
increase in the region Nf = 6 — 8, while it is approximately constant in the region Nf < 4. We have 
discussed this qualitative change for Nf > 6 and a possible relation with a preconformal phase. The 
ratio Tc/Aief is a decreasing function of Nf. This behaviour is consistent with the result obtained in 
the functional renormalization group analysis ||5|]. Next steps of the current project involve a scale 
setting at zero temperature by measuring a common UV observable. 



6 



Thermodynamic Study for Conformal Phase in Large Nf Gauge Theory 



Kohtaroh Miura 



Acknowledgements 

We thank Holger Gies, Jens Braun, Michael Muller-Preussker, Marc Wagner, Biagio Lucini, 
Volodya Miransky, Albert Deuzeman, and Tiago Nunes da Silva for fruitful discussions. This 



work was in part based on the MILC Collaboration's public lattice gauge theory code []18|]. The 
numerical calculations were carried out on the IBM-SP6 at CINECA, Italian-Grid-Infrastructures 
in Italy, and the Hitachi SR- 16000 at YITP, Kyoto University in Japan. 

References 

[1] K. Miura, M. P. Lombardo and E. Pallante, arXiv:! 1 10.3152 [hep-lat]. 

[2] W. E. Caswell, Phys. Rev. Lett. 33 (1974) 244; T. Banks and A. Zaks, Nucl. Phys. B 196 (1982) 189. 
[3] V. A. Miransky and K. Yamawaki, Phys. Rev D 55 (1997) 5051 [Erratum-ibid. D 56 (1997) 3768]. 
[4] See Ref. [0], where the analytic results are summarized with references. 

[5] J. Braun, C. S. Fisher, H. Gies, Phys. Rev D84 (201 1) 034045; J. Braun and H. Gies, JHEP 1005 
(2010) 060; 0606 (2006) 024. 

[6] M. Velkovsky and E. V. Shuryak, Phys. Lett. B 437 (1998) 398. 

[7] For recent reviews, see L. Del Debbio, PoS LATTICE2010 (2010) 004; E. Pallante, PoS 
LATTICE2009 (2009) 015. 

[8] T. Appelquist, G. T. Fleming, M. F Lin, E. T. Neil, D. A. Schaich, Phys. Rev D84, 054501 (201 1). 

[9] T. Appelquist, G. T. Fleming and E. T Neil, Phys. Rev. D 79 (2009) 076010; 
Phys. Rev. Lett. 100 (2008) 171607 [Erratum-ibid. 102 (2009) 149902]. 
[10] A. Deuzeman, M. P Lombardo and E. Pallante, Phys. Lett. B 670 (2008) 41. 
[11] A. Deuzeman, M. P Lombardo and E. Pallante, Phys. Rev D 82 (2010) 074503. 
[12] A. Hasenfratz, Phys. Rev. D 82 (2010) 014506. 
[13] A. Hasenfratz, Phys. Rev. D 80 (2009) 034505. 

[14] Z. Fodor, K. Holland, J. Kuti, D. Nogradi and C. Schroeder, Phys. Lett. B703 (201 1) 348-358. 

[15] Z. Fodor, K. Holland, J. Kuti, D. Nogradi and C. Schroeder, Phys. Lett. B 681 (2009) 353. 

[16] See reviews in Ref. [Q], and references are therein. 

[17] T. Appelquist et al, Phys. Rev. Lett. 104 (2010) 071601. 

[18] MILC Collaboration, http://www.physics.indiana.edur sg/milc.html 

[19] C. Bernard et al, Phys. Rev. D 75 (2007) 094505. 

[20] M. Luscher and P Weisz, Phys. Lett. B 158 (1985) 250; Commun. Math. Phys. 97 (1985) 59 
[Erratum-ibid. 98 (1985) 433]. 

[21] M. A. Clark, PoS LAT2006 (2006) 004. 

[22] G. P Lepage, P B. Mackenzie, Phys. Rev. D 48 (1993) 2250. 

[23] S. Gupta, Phys. Rev. D 64 (2001) 034507. 

[24] H. Kawai, R. Nakayama and K. Seo, Nucl. Phys. B 189 (1981) 40. 
[25] L Braun, Phys. Rev. D81 (2010)016008. 



7 



